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The density of states of a particle in a 2D area is independent both of the energy and form of the area only at 
the region of large values of energy. If energy is small, the density of states in the rectangular potential well 
essentially depends on the form of the area. If the bottom of the potential well has a potential relief, it can 
define the small eigenvalues as the discrete levels. In this case, dimensions and form of the area would not 
have any importance. If the conservation of zero value of the angular momentum is taken into account, the 
effective one-particle Hamiltonian for the 2D electron gas in the magnetic field in the circle is the Hamiltonian 
with the parabolic potential and the reflecting bounds. It is supposed that in the square, the Hamiltonian has 
the same view. The 2D density of states in the square can be computed as the convolution of 1 D densities. The 
density of one-particle states for 2D electron gas in the magnetic field is obtained. It consists of three regions. 
There is a discrete spectrum at the smallest energy. In the intervening region the density of states is the sum 
of the piecewise continuous function and the density of the discrete spectrum. At great energies, the density 
of states is a continuous function. The Fermi energy dependence on the magnetic field is obtained when the 
field is small and the Fermi energy is located in the region of continuous spectrum. The Fermi energy has the 
oscillating correction and in the average it increases proportionally to the square of the magnetic induction. 
Total energy of electron gas in magnetic field also oscillates and increases when the magnetic field increases 
monotonously. 
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1. Introduction 

The spectrum of one-particle Hamiltonian in a bounded area is discrete. When values of energy are 
more than maximum of potential energy (in the fourth section this criterion will be improved), the dis- 
tances between levels are of the order H 2 /2mW 2ID . Here, mis the mass of the particle, Wis the volume of 
the area, and D is the dimensionality of space. When the volume of the area is macroscopic, the spectrum 
may be considered as quasicontinuous.Then, the density of states 3lo(.E) may be introduced: 

N D {E + 8E)-N D {E) 

9Jd W = — ^ — ■ (1) 

oE 

Here, Nu {E) is the number of states, energy eigenvalues of which are less than E, 8E is a small interval 
of energy which is larger than the distances between discrete levels. Function Nd[E) was considered 
mathematically rigorously in the monograph |1]. The Schrodinger equation was multiplied by 2m/h 2 , 
and the eigenvalue £ = 2mE/h 2 has the dimension of physical quantity [L~ 2 ] (for brevity let us also refer 
to it as "energy")- It was shown that the asymptotical formulae at £ — ► oo for functions No{£) are: 

S — V 

N 2 {£)< £ + 0C\/£, N 3 {£)< T £ 3/2 + 0C£, (2) 
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when the wave functions are equal to zero at the boundaries of the area. Here, \9\ < 1 and c is indepen- 
dent of e. It is important that the first terms of these formulae are independent of the form of the area. 
The functions Nd(e) take only integer values. If this is neglected when energy is great, the formula |Q} 
can be considered as 9io = dNo/de. Then, we can obtain from the formulae Q: 

S V _ 

9?2 = — , 9l3=-7^ (3) 
4n An 1 

if only the first terms of the asymptotical expansions are taken into account. Evidently, extrapolation of 
the formulae (3} to small values of energy is incorrect even when there is no potential energy. 

In the second section of this work, the density of states at small energy is considered in the absence 
of a potential energy. 

In the third section of this work the 2D density of states is considered in the potential well that is 
created by the harmonic potential and the reflecting boundaries. As is shown in the work 0], this po- 
tential takes place in the effective one-particle Hamiltonian of the electron gas in the magnetic field. The 
dependencies on the magnetic field are investigated for the Fermi energy and for the total energy of the 
gas. 

The fourth section is a mathematical supplement. The problem on the linear harmonic oscillator with 
condition zeroes of wave function at the segment ends is considered. 



2. The density of states in a square 

The derivation of the formula (2) for AT 2 (e) starts from consideration of the square with S - L 2 

J~~L 

(see monograph |1]). The states are determined by two integer numbers. The energy of state \kl) is 
£ ki - Ji 2 {k 2 + l 2 )IL 2 . Then, A^fe) will be equal to the number of junctions of the net of squares that are 
parallel to the coordinate axes and have the sides equal to unit, which fit into the interior of the positive 
quadrant of the circle with radius {Lin) y/e. This quantity differs from the area of the quadrant SelAn by 
the sum of areas of partial squares that are crossed by the circle. The second term in the formula (|2) is 
the approximate estimate of this amendment. The relative magnitude of this amendment will be smaller, 
when the quadrant radius is larger. 

Let us consider the other method of calculating the state density in the quadratic area that can be 
used for small energy values too. The 2D Schrodinger equation for a free particle, provided that the wave 
function is equal zero at the boundaries of the square, can be changed by two identical ID equations. An 
eigenvalue of the 2D equation is the sum of eigenvalues of ID equations. Therefore, let us consider the 
state density for the ID equation. 

The number of states, whose eigenvalues are less than e for ID equation, is as follows: N\ (e) = 
\L s/e/n] . Here, \a\ denotes the integer part of the number a. The density of states that is determined by 
formula is an interval function rather than a point function. The magnitude of the interval cannot be 
taken arbitrarily small. In ID, this magnitude is limited by the demand that one eigenvalue should be in 
the interval at the greatest energy £ m . Then, 

L / , -\ 2n 2n , 

-[Vem + 8- Ve^j>l, 8>[j) + — s/£^~j-V^- (4) 

In most cases, the density of states is used in integral formulae. Then, an interval function can be changed 
by a piecewise continuous stepped function or a continuous differentiable function that is determined by 
any method of interpolation. It is apparent that the consideration of peculiarities of the state density is 
meaningless. 

Let us determine the state density for the square at the values of energy that are smaller than c m . The 
interval 8 is determined by the formula {4), an d it is accepted as the unity of energy. Non-dimensional 
(e/8) eigenvalues of energy for ID equation are denoted as Ai and A2, and non-dimensional eigenvalue of 
energy for 2D equation is denoted as \i, \i — X\ +A2. Let us consider the intervals [Ai, Ai + 1] and [A2-I, A2] 
where A2 = \i- X\. The eigenfunctions of the 2D equation that are products of the eigenfunctions ID 
equations, which are related to these intervals, have the eigenvalues that are located at interval l,jtt+ 
1]. The number of these states is denoted as MzijJ- 1). In a similar way, the eigenvalues of the 2D 
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equation that are located at the interval [^,/u+2] are obtained when A2 = jit- Ai + 1. Then, the number of 
the eigenvalues of the 2D equation at the interval [[i, \i + 1] (when /i 3= 2) is: 

M 2 (jU,;u+l) = - £ Mi(A,A+l)Mi(/z-A-l,/z-A) 



' 1=0 

1 A ' 

+ - £ Mi(A,A+UMi(jU-A,jLt-A+l). 

2 A=0 



It follows from the formula lTJ) that: 



Mi(A,A+l) = ■^^(v / A+l- s/lj. 



Then, M 2 (^n+ 1) = {L 2 8 12n 2 )l{^), where 



M-i, 



p-i 



1=0 

This function of integer argument can be written as: 

I{H) = SCju + 1) - S(/i) - S(/i - 1) + Sin - 2), 

where 



S{x) = £ \/{n+l){x-n), xTsl; S(0) = 0. 

The function Six) can be computed using Euler-Maclaurin method: 

I x+ 1 \ 2 x-l _ x- 1 
S(x) = arcsin i-xyx- 



x+ 1 



12v^c 



Based on the approximate formula 

X-l 7T 



arcsin 



x+1 2 v / x+T 3(x+l) 3/2 4(x+l) 5/2 



the asymptotical formula 



n 3 
EQi) « - + ■ 



(5) 



(6) 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



2 4^ 

can be obtained. The function Z{[i) can be calculated numerically. The results are obtained from the 
formula and from the asymptotical formula l lT2t tabulated in table 1. The values of the function lifi) 
approach nl2 = 1.57096 from above. 



Table 1. The function £Qu) that is numerically calculated [formula (7)], and calculated using an asymptot- 
ical formula (12). 



A 4 


1 


10 


100 


1000 


10,000 


100,000 


1,000,000 


formula (7) 


1.81 


1.577 


1.571 


1.5708 


1.570797 


1.570796 




formula (12) 






1.65 


1.59 


1.578 


1.573 


1.5715 



For sufficiently large value \i the density of states in a square is described by formula: 

„ M 2 {u,u+l) I 2 L 2 ?,L 2 S 
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The first term in this formula coincides with the common expression for 2D system. It can be obtained by 
differentiation A^Ce) [formula (2)], where S-L 2 . This formula can be used for derivation of the density of 
states for a flat geometrical figure of arbitrary shape (see monograph I1D . In this process, L 2 is changed 
by the figure area S and amendments are proportional to £~ 1/2 , i.e., they alter the second term in the 
formula (l3i . Therefore, the second term depends on the figure shape and on the interval magnitude S. 

The sign of the second term also depends on the figure shape. In the square this term is positive, 
i.e., when the energy increases, the density of states decreases tending to the constant value from above. 
This is explained by the fact that in ID, the state density increases when the energy decreases. The eigen- 
functions in a circle are the Bessel functions of the first kind J n {r s/e). The eigenvalues in this case are 
E nk = j 2 k IR 2 where R is the circle radius, and is the null of the function /„ that has the number k in 
the order of increasing. There is no formula that describes these nulls when their numbers k are small, 
but it is known that the distances between nulls increase when their numbers decrease. Therefore, the 
density of states should decrease when the energy decreases, and the amendment should be negative. 

In fact the spectrum at small energy values is formed by the potential relief of the bottom of the poten- 
tial well. The distances between energy levels are determined by parameters of this relief. Therefore, this 
spectrum cannot be considered as quasicontinuous. The density of states in this case can be described by 
the set of 5-functions. By virtue of the fact that the determining factor is the potential relief, it is believed 
that the figure shape is of no significance. Then, it may be helpful to obtain the state density for the square 
in the whole region of energy values. 

If the 2D Schrbdinger equation with the potential energy can be solved by separating the variables in 
the Cartesian coordinates, then the 2D density of states can be obtained. Every interval [£i , £1 + de] on the 
axis of energy of ID states contains 9li(£i)d£ states, whose wave functions are y/(xi). Products of these 
functions with the wave functions y/(x2i that have the energy values in the interval [e2,£2 + de], where 
£2 = £ - Eu are the wave functions of the 2D states, the energies of which are in the interval [e,e + df]. 
The number of these states connected with the energy value £1 is: 

dM 2 (£, £1) = 9li (£i)9ti (£ - £i)d£d£. (14) 

Then, the density of states in the square is: 

9t 2 (£) = dM2t£) = f$Ki(a)$Ri(£-a)da, (15) 
d£ J 


i.e., it is a convolution of ID densities of states. 

The form of the spectrum in the region of small energy values can play a significant role depending 
on the form of the potential relief of quantities that are determined by integral formulae. 



3. The density of states and energy of 2D electron gas in the magnetic 
field 

It is shown in the work |2] that the statistical operator of the electron gas in the magnetic field is de- 
fined by effective Hamiltonian that is the sum of the same one-particle Hamiltonians. Each one-particle 
Hamiltonian describes the particle in the potential well with a harmonic potential and reflecting bound- 
aries. Electrons interact with each other and with the neutralizing background. The electron density in 
the magnetic field should be distributed in such a way as to shield the harmonic potential. It is shown 
in the work 1 2] that this shielding in the circle with radius R leads to renormalization of the electron 
charge e r - e \/oqIR, where (— e) is the electron charge, ao is the Bohr radius. The residual harmonic po- 
tential is proportional to to 2 where a> - e r H/m is the cyclotron frequency, H is the magnetic induction, 
m is the electron mass. 

Let us suppose that in the square with the side 2L and zero of coordinate system in the center, the 
effective one-particle Hamiltonian with the symmetrical gauge also has the view: 

H 2 ( d 2 d 2 \ ma) 2 , 9 

h = h-7 + + (* + y) ■ (16) 

2m \dx 2 dy 2 8 v y ' 
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Here, e r = e V ' a^lh. Separating the variables and multiplying by 2m/7z , we obtain two identical ID equa- 
tions: 



„ 2mE v m 2 d) 2 xl 
^+——^-——^ = 0, v =l,2. 



H 2 r Ah 2 

The boundary conditions are: 

ifr(±L) = 

The problem on a linear oscillator with the boundary condition ( fl&t is considered in the fourth section. 
In this case, the ID density of states is as follows: 



(17) 
(18) 



SRi(e) = X>(£-£„) + 0(£-£ b ) — -, 

n=0 71 V £ 



n ■■ 



4mcoL 2 
n 2 h 



A = 



h 



(19) 



Here, the first term is the spectrum of the linear oscillator without reflecting boundaries. (We change n + 
j{n) to n because j{n) are practically everywhere small). For £ > £ b , the density of states is the density of 
quasicontinuous spectrum. Boundary magnitude n — hq is determined approximately, but it is significant 
that no is proportional to the magnetic induction and it is an integer. 

The density of states of Hamiltonian fl6) can be obtained by the formula fl5l It consists of three 
regions. At the smallest energy, the spectrum is discrete: 



9?2d(£) = oUb+--£ £ n8{e-An). 

^ / n=l 



(20a) 



In the intervening region, the density of states is the sum of a piecewise continuous function and the 
density of a discrete spectrum: 



9?2 P c(£) 



0[£---£ b j0(2£ b -£) 

2no+l 2L " E 

x < {2n + 2-n)8[E-An) + — 



n=na+2 



7T 



(=0 



£ - A 2 + - 

2 



-1/2 



n e = ["— j-(« + l)- 
In the region of great energy values, the density of states is a continuous function: 



f 2L n ° 

9t 2c (£) = 0(£-2£ b W — £ 



A i + 



-1/2 



2L Z / 2£ b 
H =- arcsin 1 



(20b) 



(20c) 



These formulae are illustrated in figure 1. Let us represent the states by the points at the first quadrant. 
The Cartesian coordinates of the point are the eigenvalues of ID Hamiltonians that are the components of 
the 2D Hamiltonian. The density of states in the point is the product of ID densities in the projections of 
the point multiplied by d£ [see the derivation of formula (15)]. All the states with the same energy £ - a 
are located at the line segment that cuts off the intercepts equal to a on the axes. The 2D density of states 
that have the energy a is the sum of the densities over all points of this line segment. The coordinates 
of the states that are located in the square 0£\]jDE2b or on its sides are as follows: £; = A(« ; - + 1/2). 
Therefore, the states that have energy eigenvalue £ =S £ b + A/2 (for example the point A) form the discrete 
spectrum [formula (20a)]. The degeneracy multiplicity of the discrete state with energy £ n = nA, {n = 
1, 2, ... no + 1) is equal to n. These degenerate levels are transformed in zonule, if the quantities j{ri) 
are taken into account, but this broadening is small everywhere except the immediate neighborhood 
of £ b . Among the states with energy £ b + A/2 < £ < 2£ b [the section b\b2, formula (20b)] there should 
be such ones that are located outside the square 0£u,D£2b- One of the projections of such a state is 
located in the region of the continuous ID spectrum (the point B). Therefore, the formula (20b) consists 
of a discrete part and a piecewise continuous part. The degeneracy multiplicity of the discrete state with 
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Figure 1. The illustration of calculation of the density of states for 2D electron in the magnetic field in a 
square. 

energy e„ = rcA, [n = rio + 2,...2rio + 1) is equal to 2no + 2- n. When the state energy e > 2£ b [the section 
c\Cz, formula (20c)] the point C' is similar to the point B, and for the point C, the both projections are 
located in the regions of a continuous ID spectrum. 

Let us consider the weak field case, when the Fermi energy £p > 2e\,. If the density of states is deter- 
mined by averaging over the interval A, and the amendments that take the form of the area into account 
are neglected, the density of states 912(e) should be like the one obtained in the work |2]: 



5R 2 (e) = 



m 2 u) 2 S 
4nh 2 



m 2 t) 2 S 
4nh 2 



_S_ 

4~U 



(21) 



The integral of the density of states that determines the number of states whose energies are smaller than 
e, can be represented as follows: 



JV(e>2e b ) = J m 2 (a)da + J [$K 2 (a) -3t 2 (a) | da, 
o o 



(22) 



where 9i2(o;) is taken from the formulae (20a, 20b, 20c). The equation for calculation of the Fermi energy 
£p is obtained by equating AT(£p) to the total number of electrons Nq. (For the sake of simplicity, the 
spin and the Pauli paramagnetism is not considered). This equation can be considered as the implicit 
definition of the function e?(H). If the second term in the formula (22) is neglected, then: 



£ F (H) = £ F o + £fi (H) = £fo + ■ 



1 lmwL\ 2 



2n\ h 



1 7 

£fo + — {LAY 

271 



(23) 
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where epo = An No IS is the Fermi energy in the commonly used theory when there is no magnetic field. 

The Fermi energy depends on the magnetic field due to the dependence on the magnetic induction of 
two parameters: £f(H) = £f [A(H), no{H)}. These dependencies are qualitatively different: A(cu) = mw/h 
is the linear function of H that can take on any values, and no{(i>) = \AmwL 2 ln 2 h\ can be only integer. 
Therefore, when the magnetic induction varies in an interval, in which the parameter no does not vary, 
the singular points of the function 9*2 M [formulae (20a)-(20c)] move continuously and e^{H) = £p IMH)] 
varies continuously. When the variation of the magnetic induction changes the parameter no, the number 
of the singular points changes, the spectrum reconstructs, and £p varies non-continuously. Formula (23) 
that describes the function £?{H) when the density of states is smoothed, should be supplemented by the 
oscillatory term £f2(A, no). By integration 9i2(a) between the limits and £p > 2ft, we obtain: 



N(e F > 2e b ) = ni 



ALsfK 



2 -f (££-,) 



1/2 



2L 2 



2A.no \ i 

eparcsin 1 2 \J Aho(£f _ Ano) + 2A«o 



(24) 



If n (H) - n {H + h) = +1 then A(H) - A{H±h) = + {ti 2 IAL 2 ). Then, the continuous change of £ F on this 
interval is: 



d£ F / n z 

O c £v = -r-r ±777 
dA \ AL 2 



dN (dN 



dA \ de 



AL 2 



1 dN I n l 



7T 



A 



A n 



(25) 



This change is much larger than the change that is described by the formula (23). It describes the oscilla- 
tions of £f2(A, no). At the end of the interval of the continuous change of £f2(A, no), when no changes by 
+ 1, the jump of the function £p is: 



Sje F = : 



d£p 
dn 



dN IdN 
dno \c?£f 



71 dN _ 71 /72rj£FO'\ 



(26) 



When the magnetic field varies, the Fermi energy oscillates with the amplitude {HIL) \ / 2Efolm and the 
period SH = n 2 hlAL 2 e r . If L 2 - 10" 4 m 2 and e r - 10" 4 e, this period is of the order of 10" 6 T. The 
monotonous change of the Fermi energy with the magnetic field is described by the formula (23). 

It can be similarly proved that the energy of the electron gas in the magnetic field is described by 
formula: 

S 2 p e 2 H 2 S 3 e r 4 H 4 



E=E +- 



+ E n 



(27) 



167rm I53&7i 3 h 2 m 

Here, Eq is the energy of the electron gas in the absence of the magnetic field, po = No IS is the density of 
gas. This formula differs from the energy of gas that was calculated in the work |2] by an ultimate term 

-Eosc that describes the oscillation of the energy with the amplitude (4/3^) J {2mL 2 lh 2 )E^ Q . 

We take into account that the number of discrete levels no is an integer and obtain the characteristics 
of the gas that have the oscillating dependence on the magnetic field. In the commonly used theory, the 
degenerate multiplicity of equidistant levels d{H) is an integer and is the same for all levels. This quantity 
differs from no only by the numerical coefficient that is of the order of unity. When the magnetic field, 
for example, increases, the gas energy should increase linearly until d is constant. When d increases by 
1, the \Nold~\ electrons should drop from the top level and the energy decreases by jumps. If the number 
of electrons on the top level is less than \No I d] , the jump amplitude should be smaller. These decreased 
jumps should recur, and only this oscillation is considered in the common theory. The 2D electron gas 
in the magnetic field has been considered in the monographs |3, 4]. The fact that d{H) is an integer 
is not taken into account in these works. In the work [2] it was shown that the system of equidistant 
degenerated levels cannot be a correct description of the one-particle spectrum of the electron gas in the 
magnetic field because in this theory, the angular momentum conservation and the Coulomb interaction 
are not taken into account. 
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4. The linear oscillator with zero boundary condition 

Let us change the variable in the equation (17) z — x vA and designate u = -2Elha>. Then, the equa- 
tion obtains the form of a standard equation for the function of parabolic cylinder (see handbook |5]): 



* 2 tfr i - 2 
dz 2 \ 4 

The even and odd solutions of this equation are: 



dV (z 2 

l — + u\y/ = 0. (28) 



, z 2 \ (u 1 1 z 2 

z 2 \ I u 3 3 z 2 



Vo = A ze W {--^- + (29) 

Here, <&[a, c; t) is the degenerate hypergeometric function (DHF), A e ( ] are the normalization constants. 
The middle of the line segment is zero of the coordinate. The length of the segment is 2L. To satisfy the 
boundary condition (18), the eigenvalue u should be such that 



,„ 11,) (u 3 3 \ , mwL 2 

<P - + -,-;z 2 =0, or O - + -,-;z Uo; z 2 = . (30) 

'2 4 2 1 \2 4 2 L L 2H 



This problem could not be consequently considered because in the description of nulls of DHF in all 
mathematical handbooks (see for example 1 5],| 6]) an inaccuracy takes place. It is proved that, if a < and 
c> 0, the number of nulls of the function <£>{a, c; t) is equal to (- a), if it is integer, and is equal to \— a] + 1, 
if {-a) is non-integer. It is also proved that nulls are described approximately by the formula: 

W * C) = 2^4^" (31) 

if \a\ » 1. Here, £,i{a,c) is the null of DHF that has the number i < \-d\ + 1 in the order of increasing, 
] 2 _ x . is the square of the respective null of the Bessel function of the first kind J c -\ [x). If this were so, the 
greatest nulls of the functions that are the solutions of the considered problem should be linear functions 
of the eigenvalue u. Then, the smallest eigenvalues should be ~ mtoL 2 12H, i.e., the spectrum should begin 
from very large values of energy. In the work 1 7] (see also |2]) it was shown that the formula (31) describes 
each null of the DHF only when [-a) is integer. Then, the Kummer power series that describes DHF 
terminates at the term with number I - a, and the DHF in the formulae (30) are proportional to the 
Laguerre polynomials L^^ 2 {z 2 ). The number of nulls of these polynomials is k = -a. All these nulls are 
real, positive, simple and are described by the formula (31). Let us consider the DHF <£(a, 1/2; t) when 
a — —k — j, where k is an integer and < y < 1 . Then, this DHF has k + 1 nulls of which k nulls come 
out of the nulls of Laguerre polynomial that are diminished by quantities which are proportional to y. 
These nulls are described by the formula (31). The null that has the number A; + 1 is the largest, and its 
dependence on y should possess the following properties: 

lim[f fc+1 (-fc-y,l/2)] -oo, lim [f fc+1 (-fc-y, 1/2)] =f fc+ i(-ifc- 1,1/2). (32) 

y— *0 j— >1 

To calculate k-j, 1/2) when y is small, the DHF <5?[-k— y, 1/2; t) should be changed by the asymp- 

totic expression. The Kummer power series would be represented as follows: 

0> (- k - y, 1 12; t ) = Pi 12 (y, t ) - y 12 (y, t) - y (1 - y) T\ 12 (y, t ) . (33) 

Here, P^. /2 (y, t) is a polynomial that can be obtained from the Laguerre polynomial L~ 1/2 (f) by changing 
k — ► k + y in its coefficients, yQ^ /2 (y, t) is the next term of the Kummer series that is proportional to y, 
y (1 - y) T^ l2 {y, t) is the remaining infinite series that also is proportional to y. It can be shown that when 
t is large and y is small 

y(l-y)r fc 1/2 (y,f)=y(-l) fc (l + y) S; v^(f)""" 1/2 exp(f). (34) 
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Here, [b] ,■ = b(b + 1) ■ • • (b + i - 1) is the Pochhammer symbol. The second term in the formula (33) can be 

fc 



neglected. The polynomial Pi 12 (y, t) can be changed by the last term, when t is large: 



* {H2) k k\ (l/2) fc fc! 

We set the obtained approximate expression for DHF equal to zero. This equation defines the quantity y 
as the function of ^ k+ i{-k-y,\l2) = to: 

7i/ 2 (fc,<b)= fc!r ( fc+ ^) ?o fc+1/2ex P(- ? o)- (36) 

For the boundary condition to be satisfied by the largest null of DHF 3>(-fc-y,l/2; t), this null should be 
larger than f fc +1 (-fc- 1, 1/2), i.e., 

ma>L 2 1 , 
*b = > r,,,!,,, (37) 

2/i i+4(fc+ir- 1/2 ' fc+1 

7-1/2,/ - f(i _ 1/2), where Z is integer. Therefore, the boundary condition can be satisfied by the largest 
null of DHF, if the value of k is not greater than ko: 



fco : 



2ma>L 21 
7i 2 H 



(38) 



This is the approximate formula, but the operation of taking an integer part emphasizes that ko is an 
integer, and when the frequency changes, this quantity does not change continuously and takes only 
integer values. The eigenvalues for the even solutions: 

E e = hwUc e + i+yi/ 2 (fce)J. (39) 

If k > ko, the boundary condition can be satisfied by one of the nulls of DHF that is described by the 
formula (31) when {-a) > ko- To calculate these values a, i.e., the eigenvalues of energy, we use the 
first term of the expansion DHF over the Bessel functions (see monograph |6]). This expansion is rapidly 
convergent, when \a\ is large. We obtain the following even solution of the equation (28): 

1 1\ I uz 2 \ 114 x 

y e ttA e Y - --J" 7-1/2(2 V-wz 2 /^ = A;Cos(z ^2E e lhti>^ =4 e cos^ s/2mE^. (40) 

The eigenvalues of the energy for which these wave functions are equal to zero at the ends of a segment 
are obtained: 

n 2 H 2 I l\ 2 7i 2 H 2 

£ e = ^Ue + - = — — r(2fc e +l) 2 , k e > k . (41) 



2mL 2 { 2 J 8mL 2 
A similar computation for the odd wave function leads to the results: 

3 



73/2 (Mo) = 



l,W\k+ ^ 



f 2fc+3/2 exp(-£ ), £o=Mfco + 7+73/2(fco)h (42) 



4 



- s/2mE \, E =- — —k 2 = ——-{2k ) 2 , k o >k . 

The formulae for the spectrum can be unified, if it is taken into account that the amendments y can be 
neglected everywhere except the immediate neighborhood of ko: 

hcot l\ n 2 H 2 2 

E„ = — \n+-\, n<n = 2k ; E n = - — -rn , n>n . (43) 
2 \ 2} 8mL z 

If in the equation (17) we change a> = 2a>o, this equation will gain the form of a standard equation for 
the linear oscillator with frequency cjq- If L — «■ 00, the boundary condition (18) should be changed by 
the requirement of normability of the wave functions. Then, the obtained solution of the problem turns 
into the common solution for the linear oscillator. When n > no, the approximate solutions (40) and (42) 
are the common solutions for the particle in the rectangular potential well. However, the energy E^ = 
hti)rio/2 = 2ma) 2 L 2 / n 2 that is the boundary between the two kinds of solutions does not coincide with the 
value of the potential energy at the boundary of the area as might be expected from the quasiclassical 
consideration. 
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rycTi/ma oflHonacTUHKOBi/ix CTamB /yin 2D e/ieicrpoHHoro ra3y y 
MamiTHOMy no/ii 

I.M. ,4y6poBCbKi/iM 

iHCTHTyT MeTa/iocf>i3HKi<i, 6y/ibB. BepHaflCbKoro 36, KmTb 03680, YKpaiHa 

TycTWHa CTaHiB nacTUHKH y 2D 06/iacri He 3a/ie>KHTb Bifl eHeprii i c|>opMW 06/iacTi Ti/ibKW npn Be/iwKnx 3Ha4eH- 
hax eHeprii. npn Ma/iiffi eHeprii rycn/iHa CTaHiB y npnMOKyTHifi noTeHLjia/ibHii/i HMi cyTTEBO 3a/ie>KHTb Bifl (fiopMH 
o6/iacri. Akluo aho noTeHLjia/ibHOi amm Mae noTeHLiia/ibHufi pe/ibecf>, to BiH Moxe BH3HanaTH Ma/ii B/iacHi 3Ha- 
neHHfl eHeprii hk flucKpeTHi piBHi. y L(bOMy BunaflKy po3Mipn i (fiopMa 06/iacri He MaioTb 3HaneHH?i. Akluo 
npuwiviaTn flo yBarn 36epexeHHH Hy/ibOBOro 3HaneHH?i KyTOBOro MOMeHTy, e^eKTWBHwPi oflHOHacTUHKOBnii 
TaMi^bTOHiaH fl/ifl 2D e/ieiapoHHOro ra3y y MarHiTHOMy no/ii y KO/ii e raMi/ibTOHiaHOM 3 napa6o/ii4HWM noTeH- 
L(ia.n0M i Bifl6wBai0HWMH rpaHWL4«MW. npwnycKaeTbcn, luo y KBaflpaTi TaMi/ibTOHiaH Mae TaKwfi caMnti bh™/). 
2D rycTHHa CTaHiB y KBa,qpaTi Moxe 6y™ o64wc/ieHa hk 3ropTKa 1 D rycTHH. 064naieHO rycTHHy CTaHiB 2D e/ie- 
KTpoHHOro ra3y y MarHiTHOMy no/ii. BoHa CK/iaflaeTbcn 3 Tpbox 06/iacreffi. Ko/iw eHeprii Ma/ii, cneKTp e flucKpe- 
thum. y npoiviixHiPi 06/iani rycTHHa CTaHiB e cyMOio npoMixKOBO-HenepepBHOi cpyHKi4ii' i rycTWHW flwcKpeTHOro 
cneKTpy. npn Be/11/iKwx 3Ha4eHHnx eHeprii rycTWHa CTaHiB e HenepepBHOio tj)yHKL(ieK) eHeprii. OflepxaHO 3a- 
ziexHiCTb eHeprii <t>epMi Bifl MarHiTHOro no/in, ko/ih none e aia6Ki/iM i eHeprm <t>epMi 3HaxoflWTb«i b 06/iar.Ti 
HenepepBHOro cneKTpy. EHeprii cpepMi Mae ,qo,qaHOK, nv.v\v\ ocLju/iroe i, b cepefl+HbOMy, 3pocTae nponopL(ii7iHO 
KBaflpary MamiTHOi iHflyKi^ii. riOBHa eHeprii e/ieKTpoHHOro ra3y y MarHiTHOMy no/ii TaKOx ocuy\i\K>e. i 3pocrae, 
ko/ih MarHiTHe none mohotohho 36i/ibLuyeTbCfl. 

K/iioiOBi cnoBa: rycmHa CTaHiB, ejieKTpoHHHfi ra3, MarHiTHe none, enepreivmHwiA cnexip, eHeprin 0epMi, 
noBHa eHeprin 
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